Coexistence of BCS and BEC-like pair structures in halo nuclei 
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We investigate the spatial structure of the two-neutron wave function in the Borromean nucleus 
^^Li, using a three-body model of ^Li -|-n -|- n, which includes many-body correlations stemming 
from the Pauli principle. The behavior of the neutron pair at different densities is simulated by 
calculating the two-neutron wave function at several distances between the core nucleus ®Li and the 
center of mass of the two neutrons. With this representation, a strong concentration of the neutron 
pair on the nuclear surface is for the first time quantitatively established for neutron-rich nuclei. 
That is, the neutron pair wave function in '^^Li has an oscillatory behavior at normal density, while 
it becomes a well localized single peak in the dilute density region around the nuclear surface. We 
point out that these features qualitatively correspond to the BCS and BEC-like structures of the 
pair wave function found in infinite nuclear matter. 
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Pairing correlations play a crucial role in many 
Fermion systems, such as liquid "^He, atomic nuclei, and 
ultracold atomic gases [1-3]. When the attractive inter- 
action between two Fermions is weak, the pairing cor- 
relation can be understood in terms of the well-known 
Bardeen-Cooper-Schrieffer (BCS) mechanism[l], showing 
a strong correlation in the momentum space. If the inter- 
action is sufficiently strong, on the other hand, one ex- 
pects that two fermions form a bosonic bound state with 
condensation in the ground state of a many-body system 
[4-8]. The transition from the BCS-type pairing corre- 
lations to the Bose-Einstein condensation (BEC) takes 
place continuously as a function of the strength of the 
pairing interaction. This feature is often referred to as 
the BCS-BEC crossover. 

Recently, exploiting the Feshbach resonance with 
which the strength of effective interaction can be arbi- 
trarily varied, the BCS-BEC crossover has been experi- 
mentally realized for a gas of ultracold alkali atoms [9- 
11]. This has stimulated a lot of subsequent works, not 
only in condensed matter and atomic physics [8] , but also 
in nuclear and hadron physics [12, 13] (see also Ref. [14]). 

Neutron-rich nuclei may manifest both BCS and BEC- 
like features. These nuclei are characterized by a dilute 
neutron density around the nuclear surface, and one can 
investigate the pairing correlation at several densities [15], 
ranging from the normal density in the center of the nu- 
cleus to a dilute density at the surface. In this connec- 
tion, it is worth while to mention that Matsuo recently 
investigated the spatial structure of neutron Cooper pairs 
in low-density nuclear and neutron matters, and found 
the BCS-BEC crossover behavior in the pair wave func- 
tion although the BEC limit is not reached because two 
neutrons are not bound in free space but only form a 
low-lying virtual state (see below) [12]. In Ref. [14], 
proton-neutron (T=0) Cooper pairs were also studied in 
the same context. The strong density dependence of the 



nuclcon-nucleon pseudo-potential, as well as the Pauli 
principle, are responsible for the crossover phenomenon. 

In this paper, we investigate the implication of the 
BCS-BEC crossover in finite neutron-rich nuclei. To this 
end, we particularly study the ground state wave func- 
tion of a two-neutron halo nucleus, ^^Li. This nucleus 
is known to be well described as a three-body system 
consisting of two valence neutrons and the core nucleus 
^Li [16-21]. Since both the n-n and 7i-^Li two-body sub- 
systems are not bound, the ^^Li nucleus is bound only 
as a three-body system. Such nuclei are referred to as 
Borromean nuclei, and have attracted a lot of attention 
[22, 23] . A strong di-neutron correlation as a consequence 
of the pairing interaction between the valence neutrons 
has been pointed out in -'^^Li [17, 19], which has recently 
been confirmed experimentally in the low-lying dipole 
strength distribution [24]. This di-neutron correlation has 
a responsibility for the BEC-like behaviour in infinite nu- 
clear matter, and thus, despite there is only one neutron 
pair, ^"'^Li provides optimum circumstances to investigate 
BCS and BEC-like features in finite nuclei. 

In order to study the pair wave function in ^^Li, we 
solve the following three-body Hamiltonian [18, 19], 



H = Kc{l) + hnc{2)+Vnn + 
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where m and Ac are the nucleon mass and the mass num- 
ber of the inert core nucleus, respectively. h„c is the 
single-particle Hamiltonian for a valence neutron inter- 
acting with the core. We use a Woods-Saxon potential 
with a spin-orbit force for the interaction in hnc- The 
diagonal component of the recoil kinetic energy of the 
core nucleus is included in hnc, whereas the off-diagonal 
part is taken into account in the last term in the Hamilto- 
nian (1). The interaction between the valence neutrons 
Van is taken as a delta interaction whose strength de- 
pends on the density of the core nucleus. This kind of 



2 



pseudo-potential has been standard for nuclear pairing, 
see e.g., Refs. [17, 18]. Assuming that the core density 
is described by a Fermi function, the pairing interaction 
reads 



1 + exp[(i? - Rp)/ap 



(2) 

where R = |(ri +r2)/2\. The density-dependent term 
is repulsive, and the strength of the interaction becomes 
weaker for increasing density. We use the same value for 
the parameters as in Refs. [18, 19], in which Rp ~ 2.935 
fm in the density-dependent term. 

The two-particle wave function \E'(ri,r2), where the 
coordinate of a valence neutron from the core nucleus is 
denoted by r^, is obtained by diagonalizing the three- 
body Hamiltonian (1) within a large model space which 
is consistent with the nn interaction, Vnn- To this end, 
we expand the wave function ^(ri,r2) with the eigen- 
functions of the single-particle Hamiltonian hnc- In the 
expansion, we explicitly exclude those states which are 
occupied by the core nucleons, as in the original Cooper 
problem[l]. The ground state wave fimction is obtained 
as the state with the total angular momentum J — 0. 
We transform it to the the relative and center of mass 
(c.o.m.) coordinates for the valence neutrons, r = ri — r2 
and R= {ri+ r2)/2 (see Fig. 1) [25-27]. To this end, 
we use the method of Bayman and Kallio [28]. That is, 
we first decompose the wave function into the total spin 
S=0 and 5=1 components. The coordinate transforma- 
tion is then performed for the 5=0 component, which is 
relevant to the pairing correlation: 



*^="'(ri,r2) 
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hir,R) [rL(r)yi(H)](0") |x5=o) 



(3) 

where |xs=o) is the spin wave function. 

We apply this procedure to study the ground state 
wave function of the ^^Li nucleus. We first discuss the 
probability of each L component in the wave function. 
Defining the probability as 



Pl = 



r^dr 



R^dRlfUr, 



(4) 



we find Pl = 0.578 for L = 0, 0.020 for L = 2, and 0.0045 
for L — A. The 5=0 component of wave function is thus 
largely dominated by the L = configuration [16]. The 
sum of the probabilities for L = 0, 2, and 4 components 
is 0.6025, that is close to the 5 = probability in the 
total wave function, 0.606 [18, 19]. 

Figure 1 shows the square of the two-particle wave 
function for the L — component. It is weighted with a 
factor of r^R^. One can clearly recognize the two peaked 
structure in the plot, one peak at (r, R) — (2.2, 3.4) fm 
and the other at {r,R) = (4.4,1.8) fm. These peaks 
correspond to the di-neutron and the cigar-like configu- 
rations discussed in Refs. [17, 19, 22], respectively. Notice 
that the first peak is located at a small relative distance 
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FIG. 1: (Color online) A two-dimensional (2D) plot for the 
ground state two-particle wave function, R^lfL^oir, R)f , 
for ^^Li. It is plotted as a function of the relative distance 
between two neutrons, r, and the distance between the cen- 
ter of mass of the two neutrons and the core nucleus, R, as 
denoted in the inset. 
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FIG. 2: The ground state two-particle wave functions, 
r'^R^\fL=oir, R)\'^ as a function of the relative distance be- 
tween the neutrons, r, at several distances R from the core. 
The solid lines correspond to the two-particle wave functions 
of ^^Li, while the dashed lines denote those of ^^C. Notice the 
different scales on the ordinate in the various panels. 



between the neutrons and that the corresponding config- 
uration is rather compact in the coordinate space. 

The L = wave functions of ^^Li for different values of 
R are plotted in Fig. 2 (solid line) as a function of r. The 
three-body wave function has not been presented in this 
way as far as we know, although the coordinate system 
(r, R) has been employed in several previous calculations 
[16, 22]. For comparison, those of ^^C are also shown 
by the dashed lines with arbitrary scale. Since we con- 
sider the density-dependent contact interaction, (2), this 
is effectively equivalent to probing the wave function at 
different densities. Let us first discuss the wave function 
of ^-'^Li. At i? = 0.5 fm, where the density is close to the 
normal density pQ, the two particle wave function is spa- 
tially extended and oscillates inside the nuclear interior. 
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This oscillatory behavior is typical for a Cooper pair wave 
function in the BCS approximation, and has in fact been 
found in nuclear and neutron matters at normal density 
Po (see Fig. 4 (f) in Ref. [12] as well as Fig. 4 in Ref. 
[14]). As in the infinite matter calculation[12], the two- 
particle wave function has a significant amplitude outside 
the first node at 2.4 fm. This is again a typical behavior 
of the BCS pair wave function. Notice that the core nu- 
cleus was assumed to be a point particle in Ref. [16], and 
the oscillation of the pair wave function due to the Pauli 
principle is not seen there. As R increases, the density 
p decreases. The two-particle wave function then gradu- 
ally deviates from the BCS-likc behavior. At i? = 3 fm, 
the oscillatory behavior almost disappears and the wave 
function is largely concentrated inside the first node at 
r ~ 4.5 fm. The wave function is compact in shape, indi- 
cating the strong di-neutron correlation, typical for BEC 
when many such pairs are present. At R larger than 3 
fm, the squared wave function has essentially only one 
node, and the width of the peak gradually increases as a 
function of R. This behavior is qualitatively similar to 
a local density approximation (LDA) picture of the pair 
wave function in the infinite matter [12]. 

The present results also provide a unified picture of 
the di-neutron and the cigar-like configurations in Bor- 
romean nuclei. We have seen in Fig. 1 that, for '^"'^Li, the 
former configuration corresponds to the peak around r ~ 
2.2 fm while the latter to the peak around r ~ 4.4 fm. 
These correspond to the first and the second peaks of the 
solid lines in Fig. 2, respectively (see a typical case for 
R =2.0 fm). The transition from the BCS-like behavior 
of the wave function to the BEC-like di-neutron correla- 
tion shown in Fig. 2 thus suggests that the di-neutron 
and the cigar-likc configurations arc not independent of 
each other, but rather a manifestation of a single Cooper 
pair wave function probed at various densities. 

We have confirmed, using the same three-body model, 
that this scenario also holds for another Borromean nu- 
cleus ^He as well as for the non-Borromean neutron-rich 
nuclei ^^C and ^''O. See the dashed line in Fig. 2 for 
-•^^C. The similarity with ^^Li is striking. Namely, the 
oscillatory behavior is seen at small R < 3.0 fm, while a 
single compact peak appears at i? 4.0 fm. The surface 
condensation of the Cooper pair in several neutron-rich 
nuclei has been discussed also in Refs. [15, 20], although 
these references use a coordinate system which does not 
remove the center of mass motion of two neutrons and the 
surface condensation is less evident. We should mention 
that a similar, but less pronounced, space correlation has 
already been mentioned earlier in Refs. [25, 26] for stable 
heavy nuclei. All of this indicates that the positioning of 
strongly coupled Cooper pairs with maximum probabil- 
ity in the nuclear surface is a quite common and general 
feature, that is enhanced significantly in the neutron-rich 
loosely-bound nuclei. 

The transition from the BCS-type pairing to the BEC- 
type di-neutron correlation can also clearly be seen in the 
root mean square (rms) distance of the two neutrons. For 
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FIG. 3: The root mean square distance rvms for the neutron 
pair defined by Eq. (5). It is plotted as a function of the 
distance R (Fig. 3(a)) and the density p/po of the core nucleus 
(Fig. 3(b)), where po is the normal density of infinite nuclear 
matter. The solid and dashed lines are for ^^Li and ^®C, 
respectively. 



a given value of R, we define the rms distance as 

r.^,{R) ^ sM:){R) = J {C '''^'J' ii"!' ^- (5) 

V Jo r^dr\fo{r,R)\^ 

We plot this quantity in Fig. 3(a) as a function of R. 

In order to compare it with the rms distance in nuclear 
matter, we relate the c.o.m. distance R with the density 
p using the same functional form p{R)/po = [l+exp((i?— 
i?p)/ap)]~^, as used in the nn interaction in Eq. (2). Fig. 
3(b) shows the rms distance thus obtained as a function 
of density p. The rms distance shows a distinct minimum 
at p ~ 0.4po (R ~ 3.2 fm) in "Li and p ~ 0.2po {R ~ 
4.2 fm) in ^^C. This indicates that the strong di-neutron 
correlations grow in the two nuclei around these densities. 
Notice that the probability to find the two-neutron pair is 
maximal around this region (see Fig. 1). The behavior of 
rms distance as a function of density p qualitatively well 
agrees with that in infinite matter (see Fig. 3 in Ref. 
[12]), although the absolute value of the rms distance 
is much smaller in finite nuclei, since they are bound 
systems. A size shrinking effect has been found also for 
a proton-neutron pair in infinite nuclear matter [29] as 
well as in an old calculation for the nucleus [27]. 

Finally, let us discuss how the di-neutron wave func- 
tion in ^"'^Li is modified when approaching the ^Li core 
from infinite distance. It is known that a two-nucleon 
system in vacuum in the ^5, T = 1 channel {L = S = 0) 
has a virtual state around zero energy. In regularizing the 
rms distance using the method of Ref. [30] , it is obtained 
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with the reahstic Nijmegcn potential [31] that the virtual 
state has an extension of around 12 fm. We therefore re- 
alize that in ^^Li, in spite of being a halo nucleus, the nn 
singlet pair shows a dramatic change from its asymptotic 
behavior. Approaching the core nucleus ^Li, it shrinks 
down to an rms distance rims of only 2.6 fm at a c.o.m. 
position of R=3.2 fm. At the same time, it has gained 
a maximum of binding. All this happens because of the 
well known Cooper pairing phenomenon. Pushing the nn 
pair further to the center, it feels the increasing density 
of the neutrons of the core with which the nn pair needs 
to be orthogonal. Therefore, approaching the center, the 
Cooper pair again looses binding and thus increases in 
size. What is surprising is that there exists such a well 
pronounced radius in the surface where the Cooper pair 
has minimum extension and highest probability of pres- 
ence (see Figs. 1-3). This seems a quite general feature 
common to many nuclei with well developed pairing cor- 
relations as shown in Fig. 2 (see also Ref. [32]). 

In summary, wc studied the two-neutron wave function 
in the Borromean nucleus ^^Li by using a three-body 
model with a density-dependent pairing force, and dis- 



cussed its relation to the Cooper pair wave function in 
infinite matter. Wc explored the spatial distribution of 
the two neutron wave function as a function of the center 
of mass distance R from the core nucleus, that allows an 
optimal representation of the physics. We found that the 
structure of the two-neutron wave function alters drasti- 
cally as R is varied, in a qualitatively similar way to that 
for the infinite matter. We also showed that the relative 
distance between the two neutrons scales consistently to 
that in the infinite matter as a function of density. These 
features show the same characteristics of coexistence of 
BCS and EEC-like behaviors found in infinite nuclear 
and neutron matters. 
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